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The temporal evolution of a perturbation of the equilibrium distribution of a 
condensed Bose gas is investigated using the kinetic equation which describes collision 
between condensate and noncondensate atoms. The dynamics is studied in the low 
momentum limit where an analytical treatment is feasible. Explicit results are given 
for the behavior at large times in different temperature regimes. 

I. INTRODUCTION 

A distinctive feature of the quantum kinetic theory of the condensed Bose gas arises from 
the correlations between the superfluid component and the normal fluid part corresponding 
to the excitations. This causes the occurrence of number-changing processes determining 
the relaxation to equilibrium when excitations collide frequently. As a consequence, in the 
hydrodynamic regime, a collision integral Cyi describing 1 2 splitting of an excitation 
into two others in the presence of the condensate is needed. The specific form of such a term 
depends on the dispersion relation u(k) for the energy of quasiparticles and on the matrix 
element M. of the effective Hamiltonian describing the interaction between them. 

In a series of papers before the experimental realization of BEC, Kirkpatrick and Dorf- 
man [l|, E| have derived a kinetic equation in a uniform Bose gas which includes these 
processes, and they have computed explicit values for the transport coefficients of the two- 
fluid hydrodynamics. More recently, Zaremba, Nikuni and Griffin [3] have extended the 
treatment to a trapped Bose gas by including Hartree-Fock corrections to the energy of the 
excitations, and they have derived coupled kinetic equations for the distribution functions 
of the normal and superfluid components. Many implications of their approach concerning 
transport coefficients and relaxation times have been thoroughly reported in Ref . 0| , mainly 
when k-^T ^> gn c , where g = Anam~ l is the interaction coupling constant in terms of the 
s-wave scattering length a, and n c is the condensate density. In this regime the dispersion 
law approaches the free particle law u(k) k 2 / (2m). In the opposite low temperature limit 
k^T <C gn, the relevant part of the excitation spectrum is u(k) m ck, where c = \J gnm~ l is 
the sound speed near zero temperature and n is the particle density. The kinetic equation 
in this regime has been derived in Ref. but its features have been less considered. Due to 
the absence of studies on the eigenvalue problem posed by the linearized kinetic equation, 
little is known about the time behavior of the solutions of kinetic equations in both regimes. 
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This paper is devoted to deriving some results on the free evolution of initial perturbations 
of thermal equilibrium as described by the three-excitation collision term near the critical 
and near zero temperature. By making appropriate approximations of the kinetic equation 
we shall see that the problem becomes analytically tractable in some limits, and hence the 
asymptotic behavior of the solutions for large time may be evaluated. In the first regime, 
k-gT S> gn c , we consider the small momentum limit where the equilibrium distribution func- 
tion obeys the condition no(k) ^> 1. By approximating the full linearized kinetic equation in 
this limit we obtain an integro-differential equation which, after integration by parts, agrees 
with linearized evolution equation of wave turbulence E Q with the appropriate indices. 
This result is not a priori obvious when one starts from the full kinetic equation, because 
of the strong singularity of the integrand. At low temperature, k-^T <C gn, we consider two 
opposite regimes no(k) ^> 1 and n (k) <C 1. By performing a low momentum approxima- 
tion, ck <C ksT, we will show that the wave turbulence framework still emerges but in a 
form not consistent with energy conservation, so an additional improvement is needed in 
order to restore this conservation law. In the opposite regime near zero temperature, where 
ck ^> k-^T, we shall consider an approximation based on the dominance of Beliaev damping 
processes. 

A notable feature of the linearized collision terms in the wave limit is the homogeneous 
dependence on the momentum. Such a property furnishes a systematic way to compare the 
relative orders of different collision integrals at low momentum. This is achieved by simply 
considering the degrees of homogeneity, since their values depend on some indices related 
to the dispersion law, the scattering amplitude and the number of bodies in the collision. 
Thus one can see the dominance of Cu over the binary collision term C*22- 

The plan of the paper is the following. In SecJTT]we discuss the general form of the lin- 
earized kinetic equation as written in terms of irreducible components of the perturbation. 
Sec JIIII is first devoted to deriving the low momentum approximation of the evolution equa- 
tion below the critical temperature. Then, once we show that the wave turbulence picture 
is recovered, we briefly review the general procedure for solving this kind of equation, we 
apply these techniques to the present context and we determine the asymptotic behavior of 
the perturbation for large time. Anisotropic perturbations are also discussed. In Sec. IIVI 
we derive the evolution equation in the thermal regime, with particular attention to the 
requirement of energy conservation, and we deal with the issue of the dominance of C\i over 
C*22. In Sec. |V]we consider the case of very low temperature and derive some asymptotic 
results specific for this regime. Section |VT] contains some concluding remarks. 



II. LINEARIZED COLLISION INTEGRAL IN THE FREDHOLM FORM 



The form of the three-excitation collision integral 

C 12 [n]= [ [ J R(fe,fe 1 ,fe 2 ) - R(k u k,k 2 ) - i?(fc 2 ,fci,fc)]d 3 M 3 fc 2 , 



(1) 
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depends essentially on the dispersion law u(k) and the matrix element M. of the three- 
excitation interaction. These quantities determine R by means of 

R{k,k x ,k 2 ) = \M{kMM)?[K^)-^i)-^{k 2 ))5{k-k x -k 2 )\ 

x[n 1 n 2 (l + n)-(l + ni)(l + n 2 )n]. (2) 

In this work we will consider two extreme regimes of the Bogoliubov dispersion law 



u{k) 



m \2m J 



1 1/2 



(3) 



near the critical temperature and near zero temperature. In the first limit k-^T ^> gn c the 
dispersion law and the squared amplitude become 

k 2 

~ 7^ + 9ric, (4) 
2m 

\M(k, kl ,k 2 )\ 2 = ^ = ^#, (5) 

where n c is the superfluid density. Since in this regime the average value of the energy of 
a quasiparticle is of order k B T the term gn c is subleading in Eq. @. In the opposite limit 
k B T <C gn we have 

u(k) « ck + ak 3 , (6) 
\ka(ui lm2 9gkk ± k 2 9ckhk 2 

\M{k,ki,k 2 )\ = = —— . 7 

The amplitudes Ai can be derived heuristically by golden rule arguments from the effective 
Hamiltonians in these regimes 0. Note that no distinction is made between n and n c 
in the scattering amplitude valid at T — > 0, because we neglect the small low temperature 
depletion. Since we are interested in the low momentum region, we shall retain the first term 
of the dispersion law. Here we will consider the homogeneous regime where the distribution 
function is independent of the position. 

The linearization of the kinetic equation dtn = Ci 2 [n] proceeds by the insertion of 

n(M) =n Q (u(k))+n Q (u(k))[l + n Q (uj(k)))x(k,t), (8) 

where no(u(k)) is the Bose distribution function without chemical potential with the energy 
given by the first term of Eqs. (j4j) or ([6]); the dimensionless function x(k,t) parametrizes 
the departure from equilibrium for the noncondensate distribution function. The linearized 
kinetic equation adopts the form 

n (l + n )^ = L[ X ](k,t) = J C(k,V) X (V,t)tPkf, (9) 

which must be completed with the linearized equation for the fluctuation of the condensate 
density 5n c , 

-6T l2 = - / L[ X ](k)— - (10) 



dt J L/XJV J (2n) 
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The kernel £ is symmetric in k and k' and, due to the Bose functions, falls off sufficiently 
as /c, k — > oo. We may then introduce the following scalar product 

(Xi,X 2 ) = J n (l + n )xi(k)x 2 (k)d 3 k, (11) 

and pose the eigenvalue problem 

J C(k, k') Xj (k')d 3 k' = -Xjn (l + n ) Xj (k). (12) 

The solution of the linearized equation would be formally written as the following series 

x(M) = E (13) 

j 

where the coefficients Cj would be determined by the initial condition. In this work we 
do not intend to follow this procedure since it is uncertain that there exists a spectrum 
of discrete eigenvalues for this problem. This doubt is based on the observation that the 
lifetime of a long-lived well-developed sound mode of energy u = ck varies continuously with 
the momentum, r(k) oc k~ 5 , where 5 = 1 at T ^ [9] or 5 = 5 at T = [10|. So, this 
property does not seem to be consistent with the discrete character of the spectrum. Instead, 
we will try to find an expression for the action of the operator L[x] at low momentum, and 
then solve the evolution equation found in this approximation. 

Formally, the operator L is decomposed in two types of contributions 

L[ x ] = -Af(k) x (k, t) + J U(k, k') X (k', t)d 3 k', (14) 

where the function U is given by 

U(k,k') = [2\M(k,k',k-k')\ 2 5(u(k)-uj(k')~uj(k-k')) 
x n(w)[l + n(w')][l + n(w - u')] + (k <-> k')} 
-2\M(k + fe', k, k')\ 2 5(uj(k) + u(k') - u(k + k')) 
x[l +n(w)][l + n(co')]n(co + co'), (15) 

and the explicit form of the coefficient M{k) will be not required. 

The conservation laws of energy and momentum imply the presence of zero modes of L 
proportional to the energy and the momentum, x(k,t) oc u(k),h(t) ■ k. In particular, 

L\pu(k)] = -M(k)pu(k) + J U(k, k')Pu(k')d 3 k' = 0. (16) 

Thus, the linearized collision operator can be written as 

nX(k,t) 



L[ X ] = JlA(k,k') 



pu(k) 



d 3 k', (17) 
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which suggests to introduce a dimensionless variable A(k,t) defined as 

a ^=My (18) 

In terms of A, the linearized kinetic equation has the form 

n (l + n )/3u{k)d t A{k, t) = J U{k, k')/3u(k') [A{k', t) - A{k, t)} d 3 k'. (19) 

The rotational invariance of the kernel U(k, k') may be exploited by expressing the pertur- 
bation as a superposition of angular momentum eigenstates 

A{k,t) = J2Am{k,t)Y lm {k). (20) 

l,m 

Using the addition theorem for spherical harmonics one arrives at the equation 

n (l + n )/3u(k)d t Ai m (k,t) = K t [Ai m ](k,t) (21) 
where the operator K\ adopts the Fredholm form 

K t [A]{k,t) = / Ui{k,k')f3uj{k') [A(k',t) - A(k,t)}k' 2 dk' 

-AnA(k, t) J™ (u (k, k') - ^f^Pj /3uj(k')k' 2 dk'. (22) 

Here Ui(k, k') is the /-coefficient in the series of Legendre polynomials of U: 

oo 

U(k,k') = J2Ui( k i k ') p i( cosd >°k>)- (23) 

1=0 

Next we derive two approximations to -K"/[^4.](A;, t) at low momentum valid near the critical 
temperature and low temperature. They must formally satisfy energy conservation which is 
achieved if 



d 
dt 



poo poo 

' n (l + n )l3u(k)A 00 (k,t)u(k)k 2 dk = / K [A 00 }(k,t)u(k)k 2 dk = 0. (24) 
o Jo 



Clearly, in view of Eq. f l22|) . the symmetry of the kernel Uo(k,k') is sufficient in order to 
accomplish this requirement, so such a property must not be lost in the approximated kernel. 

III. DYNAMICS NEAR THE CRITICAL TEMPERATURE 

A. The linearized equation at low momentum 

In this section we consider the regime k^T ^> gn c near the critical temperature. Here 
the dispersion is well approximated by u = k 2 /{2m), and the squared amplitude by \Ai\ 2 = 
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FIG. 1: The kernel Uo(k, £/)/c' 4 for kW/3/m = 1/2 near the critical temperature. For this value 
there is no appreciable difference between the graphs of the exact kernel and the approximation of 
Eq. ()25|) . The singularity at k' = k is not integrable. 

8n c a 2 m~ 2 . The low momentum approximation to Uo(k,k') is obtained by the rescaling 
k — » ek,k' — )■ ek' and the expansion around e = 0. This gives 

„ , A 128m 2 a 2 n c ( 6(k - k') 9(k' - k) \ 

The kernel Ui(k, k') has a non-integrable singularity proportional to [A; — k'l^ 1 as k' — > k (see 
Fig. [Q. In order to manage carefully this singularity, it is convenient to integrate by parts 
in the exact expression of i^[.A] given by Eq. fl22|) . It helps to write 

r-k 1 

A(k',t)-A(k,t) = dqd q A{q,t), (26) 
Jk 

and to perform the change of the order of integration. This yields 

r-k poo 

Ki[A](k,t) = / dqJ<(k,q)d q A(q,t)+ / dq J> (k, q)d q A(q,t) 

JO Jk 

-AnA(k, t) J™ fu (k, k') - ^^Pj ^(k')k' 2 dk', (27) 

where 

J<(k,q) = -9(k-q)^j [ q U l (k,k')f3oj(k')k' 2 dk', (28) 



o 

"00 



4tt r°° 

J?%q) = %-*)2|7i/ U^k'^uik^dk', (29) 
and 6(x) is the unit step function. 
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The behavior of the first term in the right hand side of Eq. (|27|) as k — > is derived by 
determining the leading behavior of J^(k,q) as both arguments k, q — > 0. With respect to 
the second term in Eq. ( 1271) it is far from obvious that the limit q — > of J^k^q) must 
be incorporated since q runs from A; to oo. The key idea emerges by noting first that the 
considered scattering amplitude involves the assumption that all three momenta are of the 
same order. It is then expected to be adequate to take the limit q — > when k — > in 
Jf"(k,q). Another more formal justification arises from the fact that Jf"(k,q) diverges as 
ln(g — k) as q — > k. Therefore, by assuming that d q A(q, t) falls off sufficiently to insure 
integrability, the main contribution to the second term in Eq. ( 1271) as k — > comes from the 
region where q is also close to zero. The leading behavior of the integral in the last term 
of Eq. ( p7|) as k — > is easily determined since the singularities of (21 + l)Uo and U\ cancel 
each other. 

The evaluation of the integrals in that limit can be carried out analytically, and they 
assume the scaling form 

j<(k, q ) ~ —nf(-)e(k- q ), 
q \qj 

j>(k, q ) ~ — n>(-)e(q-k), 
q \qj 

Explicit expressions of Hf' (x) will be given below. In this approximation the action of the 
operator Ki is written as 

Ki[A](k,t) ~ AT 2 r U x (-) d q A(q,t)— - k-%A(k,t), k 0, (32) 

Jo \qj q 

where 

U(x) = nf(x)9(l -x)+ Hf(x)6(x - 1). (33) 

Note that X Q = 0. Since the right side of Eq. ( J2T]) behaves as d t A/({3cu), the kinetic equation 
will be approximated by 

d t A(k,t) = /3u(k)Ki[A}(k,t), k^O. (34) 

The integral operator that we have obtained in Eq. ( 13 2 p is slightly different from those that 
appear in the linearization of the kinetic equation of wave turbulence. However, by assuming 
the validity of integration by parts and noting that uq ~ (/3u;(A;)) -1 at low momentum, the 
kinetic equation (T2T]) adopts the form 

d t A{k,t) = ^u{k)Ki[A\{k,t) = k- h f°°Vi (-) A(q,t)% k -> 0, (35) 

Jo \qj q 

where h = 1. The kernel Vi(x) may include a term proportional to 5(x — 1), but in our 
case this does not occur. Essentially, such an equation is the final result of the linearization 
procedure of the nonlinear equation of wave turbulence 0] , which may be obtained from the 



k, q -)■ 0, 
k, q — >■ 0. 



(30) 
(31) 
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original ([1]) by keeping the leading part of the integrand when n>l ll|. Nevertheless, in the 
absence of the large momentum part of the distribution functions assuring the fall off of the 
integrand at large momentum, the naive justification of this procedure may be inappropriate 
and a more delicate approach is required. That is the reason why we have addressed the 
linearization of the full quantum kinetic equation preserving suitable integrability properties. 

In the remainder of this section we shall focus on the solutions of Eq. ( |34|) with K\ 
approximated by Eq. (132]) . 



B. Evolution of perturbations close to T c 

In order to obtain the characteristic time scale of the kinetic equation, it is convenient 
to introduce the dimensionless momentum variable k = k^J (3/m; the scale (m//?) 1 / 2 cor- 
responds to the average value of the momentum carried by a quasiparticle in the weakly 
interacting gas. It turns out that the right side of Eq. (13"4"|) reads 



d t A(k,t) = - U 



-\d- q A{q,t)^--^A{k,t), 
q q rk 



(36) 



where r 1 sets the inverse time scale corresponding to the scattering rate involving three 
excitations in the presence of the (small) condensate. It is given by 



1 
r 



32nn c a 



\fjfrn 



(37) 



Here f3 may be replaced by the inverse critical temperature of the ideal gas (3 C = 
mC(3/2) 2,/3 /(27rn 2 / 3 ) since the above expression is linear in the small depletion n c /n. We 
set t = tjr for the reduced variable time. The explicit form of H.q[x) is derived in the Ap- 
pendix n (the computations for the cases 1 — 1,2 proceed in an analogous way). The results 



for I 



0, 1, 2 are 



H (x) 
Hi{x) 

n 2 (x) 
It 

T-2 



9(1 -x)-\n 

x 

9(l-x)-\n 

x 



1 + x' 

1 - X 2 
1 + X 

1 — X 



^6(x 
9{x- 



' x V x 4 



2 1 ^ 
^ - - In 



-9(1 

+9(x-l) 
0. 

2 - In 16, 



3 1 , 

In 256, 

2 2 



x)— In \(l + x 2 )(l 
2x 



„2\21 



2x 3 2x 



4 i 2 
X + X 



x + 1 



X 



(38) 
(39) 



(40) 

(41) 
(42) 

(43) 



Hereafter in this section we will drop the lines over the reduced momentum variables. 
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The most efficient way to manage the convolution integral of the kinetic equation is to 
use the Mellin transform with respect to momentum and the Laplace transform with respect 
to time. If J c (s, A) denotes the image of A(k,t), 

/*oo /*oo 

T(s,\)= dte- XJ Aik^k^dk, (44) 
Jo Jo 

the evolution equation becomes 

\F(s, A) = W Ul {s) [-{s - l)7{s - 1, A)] - XiF{s - 1, A) + (45) 
or alternatively 

\7{s + 1, A) = W Vl (s)F(s, A) + tt(s + 1), (46) 

where 

Wv^^-sWu^s + l)-^. (47) 

Here Wu t (s) and ^f(s) are the Mellin transforms of 7ii(x) and the initial condition A(k,0), 
while Wv ; ( s ) denotes the Mellin image of the kernel Vi(x) in Eq. ( 1351) . We will assume 
that .A(fc, 0) has compact support to assure that \I/(s) is an entire function. Let us first 
concentrate in the I = case. From Eq. ( 138)) we see that the Mellin function is given by 

W Vo (s) = -2 7E - 2tP (0 - Trcot (^) , Res G (-2,4), (48) 

where 7e is the Euler constant and ip{z) = T'(z)/T(z) is the digamma function. Note the 
reflection property Wy (l + s) — Wy (l — s), in accordance with the symmetry of U (k, k'). 
The fundamental strip where Wy (s) is analytic is Res G (—2,4), and contains two simple 
zeros at s = 0, 2. It is important for what follows to introduce the winding number of the 
Mellin function Wy l (s) as the increment divided by 2ir of the argument of Wv t (s) as s moves 
from a — ioo to a + 200 along the line Re s = a. In the I = case, one can check numerically 
that the winding number n(a) of the Mellin function vanishes when a G (0, 2). 

To proceed further one must solve Eq. (T46]) for J-(s, A). (In this paragraph we will 
assume an arbitrary non-negative value for h.) This can be accomplished by following the 



method of Ref. [12J, briefly reviewed in the sequel. These authors have shown that when 
an interval of zero winding exists within the fundamental strip, (<j_,cr + ) c (a, b), then the 
equation (T4"6"j) has a unique solution. It can be constructed using a particular solution B(s) 
of the homogeneous equation 

B(s + h) = -W(s)B(s), Res G (a, b), (49) 



whose requisite properties have been described in Ref. [12J. The most important features 
of B(s) are the meromorphic property in the strip Re s G (a, b + h), and its analyticity and 
absence of zeros in the strip Res G (a_,(j + + h). The location of the poles and zeros of 
B(s) outside that interval is dictated by the set of zeros of the Mellin function within the 
fundamental strip. They determine the asymptotic properties of the solution. In particular, 
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it turns out that the pole of B(s) with maximal real part determines the asymptotics of 
A(k, t) as k — > 0. Now, by expressing the Mellin image of the perturbation as 

•F( S ,A) = B(s)/( Sj A), (50) 
the function / satisfies the inhomogeneous difference equation with constant coefficients 

A/( S ,A)+/( S -M) = §[4, Rese(a_,a+ + h) (51) 
which has an inverse Mellin image given by 
(A + k~ h )a(k, A) = Q(k) 

i r +io ° *( s ) 



/ wi k ~ Sds > ° G (°->°+ + h )- ( 52 ) 



2-ni 

The end step is to write the convolution corresponding to Eq. (1501 and to invert the Laplace 
transform. This yields the solution 

A(k,t) = J™Z (£) exp (-q- h t) Q(q)^-, (53) 

in terms of the inverse Mellin transform of the base function 

j f>(T-\-iOG 

Z(k) = — B(s)k- S ds, Res G (er_,er + + h). (54) 

^1 J a-ioo 

Notice that when h > the asymptotic expansion of the solution ( 1531) as t — > oo requires 
the determination of the leading behavior of Q(k) as k — > oo. According to Eq. f )52|) . this 
comes from the zero of B(s) with minimal real part. 

Next we turn to the I = case. Due to the presence of the zero of W / v ( s ) at s = 0, the 
function B (s) has a sequence of simple poles at s = 0, — 1, —2, . . ., with the behavior 

then 

m ~~^WY (56) 

where B(l) is an arbitrary constant. This leads to a behavior independent of k at low 
momentum 

A ^ ~ -t^TTvT /°° exp fc °' (57) 

Note that B(l) must cancel another similar factor arising from Q(q). On the other hand, 
the zero of B(s) with minimal real part is at s = 3, and it comes from the zero of Wy (s) at 
s = 2. Using B(s + 2) = W Vo (s + l)W Vo (s)B(s) when s^lwe obtain 

B(s) ~ B(l)Wv b (l)W^(2)(s - 3), s -> 3. (58) 
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This zero determines the high momentum limit of Q(k) in the form 

«*)~- B ( i)K(T)n(2) g- (59) 

so, the long time behavior of the solution is given by 

„/, -x *(s = 3) [°°„(k\ ,~,sdq - 

A ( k >V~- *mJmw ^ / * b ex PH/?)4, f->oo, (60) 



5(1)^(1)^(2) y w " v <t 

which possess the self-similarity property 

A(k,tu) = ^(j,vj. (61) 

Finally, it is easy to derive in closed form the asymptotics in the more restricted regime 
k — > 0, t — > oo either from Eq. (157)1 combined with ( 159]) . or from Eq. (|60|) combined with (|56|) . 
The result is 

2^( s = 3) 1 
^(0)^(1)^(2)? 
where 

W^)Wv {l)WvP) = ^ (tt - In 16) . (63) 
The source term 5T\2 only receives contribution from the / = component, 

d 3 k 



^(fe^- w/ 7^ui7 ml 77^ =3 > ^ ^ 0, t ^ oo, (62) 



<5r i2 = — / n (l + n )x(M) 



(2< 



3 



1 d 



4tt 5 / 2 dt 

or using the above low momentum approximation 



n (l + n )/3oo(k)Aoo(k, t)k 2 dk, (64) 



4vr 5 / 2 <ft 7 ^(Jfe) 
m 3//2 d 



- 2^/33/Vdt^^ 1 '^' (65) 

which contains the partial Mellin transform of Aoo{k,t) for s = 1. From the convolution in 
Eq. (l53~j) one finds 

F(s,t) = B(s) exp(-t/q)Q(q)q s - 1 dq. (66) 

The substitution of the large momentum behavior of Q(k) given in Eq. fl59|) yields the 
asymptotics 

^"--"(•'fiS' Res<x (67) 



Combining this with Eq. (1651) then gives the long-time behavior of the source term 



m 



3 / 2 ^(3) 1 



12 



Since the energy of the perturbation is proportional to = 3, t), we see that 

F{s = 3,1) = *(3), t->oo, (69) 

in accordance with energy conservation. 

What about the non-isotropic perturbations? The Mellin functions corresponding to the 
I = 1, 2 harmonics are 

W Vl (s) = -2 TE + -A_-V>(^) Res G (-1,3), (70) 

Wv 2 (s) = -2 7E + lnl6+ ° 



s(s-2) 



-* cot (t) + " • Rese( - 2 - 4) - (71) 

Now Wvi( s ) has a double zero at s = 1, and the winding number k does not vanish on the en- 
tire strip of analyticity; the corresponding values are 1, (—1) if Res G (1,3), (Re s G (—1, 1)). 
For I = 2, the winding number does not assume the value zero either, but now k = 1 on the 
entire interval (-2,4), and the Mellin function does not vanish in the strip of analyticity. In 
these situations of absence of an interval of zero winding, Balk and Zakharov |12j have shown 
that the initial value problem for A(k, 0) is not well posed, signaling the non-uniqueness of 
the solution when k > on the entire interval (as in the / = 2 case), or a strong instability of 
exponential growth if k ^ (as in the I = 1 case). In terms of the Mellin function, we have 
Wv 12 (0) > 0, which at least is a sufficient condition 12| for the instability of the Rayleigh- 
Jeans spectrum n oc k~ 2 under anisotropic perturbations. Thus, it seems praiseworthy to 
assign this property to the Bose-Einstein equilibrium solution of the kinetic equation with 
a collision term given by Eqs. (j3J) and (JSJ). 



Finally, it may be instructive to compare these findings with those of Ref. [13J. The 
application of the methods of wave turbulence theory to the study of the binary collision 
term C22, approximated when uq ^> 1, showed the Mellin function W / 22(s) to be symmetric 
with respect to s = 3/2, reflecting the symmetry of the corresponding kernel Uq. At the 
same time, it was shown that W 2 2(s = 0) > while W 2 2(s = 1) = 0. As a consequence, such 
approximated form of C22 destroyed energy conservation but was consistent with particle 
number conservation. By simply shifting the minimum of W 2 2{s) one can get = 0) = 

and W22{s = 1) 7^ 0. So that another approximated form of C22, which conserves the energy 
but destroys the conservation of the particle number, appears as possible. However, given 
the reflection property of ^22(5) and the location of its zeros, it does not seem possible to 
get a low momentum approximation maintaining both conservation laws. 

IV. EVOLUTION OF PERTURBATIONS IN THE THERMAL REGIME f3ck < 1 

In this section we consider the regime u(k) k^T gn, where the energy of quasi- 
particles is uj = ck, with c = \Jgnm~ 1 . Their low-energy interactions are well described by 
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FIG. 2: The kernel Uo(k,k') in the low temperature regime when f3ck <C 1. The continuous line 
corresponds to the exact kernel for f3ck = 1/4. The dashed line is the kernel obtained with U(k, k') 
approximated by Eq. (|72p for /3ck = 1/4. 

the squared amplitude \M.\ 2 given by Eq. fl7J). By ignoring the next positive contribution 
(8m 2 c) _1 fc 3 in the dispersion law all 1A\ become equal to (21 + l)Wo, due to the collinearity of 
wave vectors in the collision integral. Therefore in this approximation, the kinetic equation 
is the same for all values of /. The low momentum limit of this kernel is 

U (k, k') ~ - 2 \ 33 ( 9 Jt^*l + B J£^$\ , Pck 0, Pch> 0. (72) 
62ir z mnp' i c i \ k k J 

A comparison between the exact and the approximated kernels is shown in Fig. El Our first 
task is to find an evolution equation consistent with energy conservation. When Eq. ( 1721) is 
substituted into Eq. (122]) . the latter turns into 

9k 3 /A 3 1 \ 
+ o ^, U-^M(M), (73) 



lirmn/3 3 c 3 \3k 3 12 j 

where we have introduced an ultraviolet cutoff A ^> k in the integral which gives the 
coefficient of A(k,t). We left for the moment this term undetermined. The well-defined 
convolution term in the above equation determines the (provisional) class of admissible 
solutions of the approximated evolution equation. Such a class would be formed by functions 
A(k,t) which grow as k — > more slowly than k~ 4 , and fall off faster than k~ 3 as k — > oo. 
In accordance with that, the integral for the Mellin transform of the gain term 

f°° 1 

- / [x^x - 1) + x- 3 6(l - x)] x^dx = — -, (74) 

Jo («-3)(s-4) 

is convergent for 3 < Res < 4. Note however that the analytic continuation of this integral 
to s = is the value of the gain term formally evaluated for A(k',t) = 1. Based on 
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energy conservation, we may insist on considering the functions A = constant within the 
admissible class of solutions, so the loss term must cancel out this value 1/12. Therefore, 
the replacement A 3 /(3/c 3 ) — 1/12— > — 1/12 produces such a cancellation. 
The approximated evolution equation obtained in this way reads 

Qk 4 f°° / h' 4 V 3 \ dV 

dMk - ') = ~ sJ^ I (f»<* " k "> + ~ *>) Aik ' '> * 

^k 4 

;A(k,t), (75) 



32irmn/3 2 c 2 ' 

which again has the form of a linearized kinetic equation of wave turbulence, where the 
degree of homogeneity is now h = —4. In fact, we will find that the above equation still 
introduces a conflict with energy conservation. This defect is repaired by adding a source 
term independent of k to the right-hand side of Eq. ( 1751) . It has the form 

q r°° 

$(t) = - — / A(k',t)k' 3 dk\ (76) 



]irmn/3 2 c 2 J 

and produces a partial cancellation of the integral term of Eq. (175|) with the result 

3fc 4 

WW = -32W/^ (M) 

9fc 4 f°° fk' 4 fc ,3 \ dk' 
+ ^— ^ / ti-u ^,*)-TT- ( 77 ) 



87rmnf3 2 c 2 J k \ k 4 k 3 J ' fc' 

This implies that the class of solutions of the improved equation must be changed. Such a 
class is formed indeed by functions which fall off faster than k~ A as k — > oo. In order to 
understand the role of the source term and how it arises, it is convenient first to discuss the 
solution of the wrong evolution based on Eq. ( I75|) . 
By introducing the function 

V(x) = -5(x - 1) - I2x- A 8(x - 1) - 12x- 3 9(l - x), (78) 

we define the Mellin function corresponding to Eq. (j75p as 

W v {s)= [ °V(x)x s ~ 1 dx = - . S{S ~ 7) AV Res e (3,4), (79) 
Jo (s-3)(s-4) 

which has zero winding in the entire strip of analyticity, and Wy{s + 7/2) = Wv(7/2 — s). 
The corresponding difference equation for the Laplace-Mellin image of A becomes 

AJ^(s-4,A) =vW v (s)J r (s,\) + ^(s-4:), Res e (3,4), (80) 

where v = 3/(327imn(3 2 c 2 ). The method used to solve Eq. (180]) is the same as before. We 
seek an appropriate solution of B(s — 4) = — Wv(s)B(s) for Res G (3,4), and we write the 
solution in the factorized form 

HsA) = B(s) J™ J^Lfc- 1 *, (81) 

1 /-Res+ioo ,w n 



2mJ Rcs _ loo B{s) 
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In view of Eq. (1801) and such a factorization, we need a solution B(s) which has neither zeros 
nor poles when Res G (—1,4). We can write the required function as 

_ r( z i £ )r(i + ^) _ 3-s sin (f ) 

[) r(4f)r(l + S) s sin (i (, + !))• 

The other solutions of the homogeneous difference equation in terms of products and frac- 
tions of T-functions, such as 

r(-f)r(i + ^) _ 3- s sin(|(i + s )) (g4) 



r (2=2) r (1 + s=Z) s sin (fa) ' 

r(¥) r H) r(i + ^)r(i + ^) 
r(¥)r(¥) r(1 + ¥)r(1 + f) 



(85) 



are not appropriate for the present problem since they have a pole at s = within the strip 
Res G (—1,4). In accordance with the factorization of J c (s, A), one finds the solution as the 
convolution 

A(k, t)= ^ Z (-) exp {-vqH) Q(q)^, (86) 
Jo \Q/ Q 

where Z(k) and Q(k) are given by the integrations in Eqs. (154]) and (I52I) with a G (—1,4). 

Thus, to derive the asymptotics of the perturbation as t — > oo, we must evaluate the low 

momentum behavior of Q(k). This is dictated by the pole of B(s)^ 1 at s — —4 which leads 

to 

Q(k) ~ _ 27/2 y- 4 ) fc 4 ; k ^ Q (g7) 

171 

If J-'lSyt) denotes the Mellin image of the solution A(k,t), the above result permits the 
calculation of the leading behavior of the quotient J r (s,t)/B(s) as t — )■ oo. It reads 



' T{s ' f) e W {-vkH)Q{k)k s ^dk 
Jo 

23/2 lt^l T (i + £ j ^ oo, Re s > -4. (88) 



B(a) 



We can see now the violation of energy conservation. In the approximation where n 3> 1 
this requires for any time that 

= 4^^ ( ' 5 = 3 - i) = - (89) 

or .F(s = 3,i) = J-"(s = 3,t = 0) = \l/(3), since \l/(3) corresponds to the energy of the initial 
perturbation. But, in accordance with the above asymptotics, the quantity J c (s = 3,t) 
assumes a value proportional to 4)i3(3)t — 7 / 4 for large time, so dE/dt ^ 0. Another 
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indication that Eq. ( !75|) is conflicting comes from the fact that, in order to check energy 
conservation, the integration of both sides after multiplication by k 2 leads to 

7 /»00 /*00 /*00 

— / A(k,t)k 2 dk = v A(q,t)q 6 dq V(x)x G dx, (90) 

Jo Jo Jo 

which does not converge for the class of functions considered up till now. 

To obtain an improved evolution equation, let us consider a factorized solution J-(s, t) 
for a different base function respecting energy conservation. If we assume that Q(k) ~ Dk a 
as k — > 0, then the limit behavior of J-"(s, t) as t — > 00 is 



F(s,t) = B(s) / expi-vk^Qi^k^dk 
Jo 

~ jB(s)T (^-^] (^)~ (S+Q)/4 , Re(s + a) > 0. (91) 

Since energy conservation requires that J-"(s = 3,t = 00) = \&(3), we see that a = —3, and 
B(s) must possess a simple zero at s = 3 in order to cancel the pole of the gamma function. 
Furthermore, to accomplish the low momentum behavior of Q(k), the prescription for Res 
in 

Q(k) = — / lMk- s ds, (92) 

2m J Rcs _ ioo B{s) 

and the analog of (l54"j) must be Res > 3. Clearly, if the base function assumes the value 

S( S ) = ^L^ Res>3, (93) 

which was discarded before, the requirement ^(s = 3,t) = \I/ (3) is accomplished for large 
time. 

Some features of the solution generated by this base function are the following. The 
inverse images of \l/(s)/£>(s) and B(s) are simply 

3 F°° 

Q(k) = A(k,0) + -J A(q,0)q 2 dq, (94) 
Z(k) = 8(k - 1) - 30(1 - k), (95) 

and the behavior of Q(k) as k — > reads 

Q(k) ~ 3^(3)£;~ 3 , k ->• 0. (96) 

Hence the corresponding solution is written as 

A(k,t) = zf^)e^>(-vqH)Q(q)^, 

= exp(-vkH)Q(k) - 3 f exp(-t;g 4 t)Q(g)^, (97) 

Jk Q 
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which for large t assumes the value 

A(k, t) ~ exp {-vk A t) Q{k), t -> oo. (98) 

Similarly, the behavior or the source term as t — > oo becomes 
1 f°° 

6Tl2 = 4^7^ J d ^ k ^ kdk 

- f exp(-vkH)k 5 Q(k)dk + 3 f ^ exp(-vqH)q 4 Q(q) f kdk 
. Jo Jo Q Jo 

POO 

/ exp(-vkH)k 5 Q(k)dk, (99) 
Jo 



v 



8vr 5 /2/3 c 

so the insertion of Eq. fl96|) produces 



^/ 4 r(7/4)vI/(3) 
We can check directly energy conservation for any time, 

OO POO 

d t A(k,t)k 2 dk = -v exp(-vk 4 t)k 6 Q(k)dk 
o Jo 

+3v I —exp(-vqH)q 4 Q(q) [ k 2 dk = Q. (101) 
Jo Q Jo 

It remains to derive the resulting evolution equation satisfied by A(k, t) given in Eq. (197]) . 
From the homogeneous difference equation evaluated for B(s) in Eq. ( 1851) it follows that 

W v {s)B(s) = -B{s - 4) = -1 + Res e (3,4). (102) 

Applying the inverse Mellin transform we obtain 

jTvg)^)^-^-!)-^*-!) 

= -A- 4 Z(fc) - 1. (103) 

The origin of the extra term proportional to k~ A is the pole of B(s) at s — 0. With this 
result at hand, the integral term of Eq. ( 175]) evaluated for A(k,t) in Eq. (]97]) becomes 

J™ V A(q,t)^- = -Ld t A(k,t)-^j\x V (~vqH)Q(q)q 3 dq 

1 a .„ ^ 3 .F(s = 4,t) 

1 12 Z" 00 

= —d t A(k,t)-- A(q,t)q 3 dq. (104) 
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The last term can be understood as a self-consistent source which restores energy conserva- 
tion. In consequence, if one uses the partial cancellation produced by this term, the corrected 
evolution equation adopts the form 

d t A(k,t) = vk 4 J™V\^)A(q,t)^, (105) 

where 

V(x) = -5(x - 1) - 12aT 3 #(l -x) + 12aT 4 0(l - x), (106) 
which gives rise to the same Mellin function as before 

w vM = - ( /X(s- 4) ' Res>4, (1 ° 7) 

but with a different strip of analyticity. The relevant interval of zero winding is now Re s G 
(7, oo) and the base function B(s) has neither zero nor poles in the corresponding strip 
Res G (3, oo). We expect that the evolution equation obtained in this way captures the 
essential features of the low momentum regime at low temperature. 

Dominance of ' C\ 2 over C 22 . To conclude this section, it is interesting to observe that 
the collision term that we have considered in this and the previous section is the prevailing 
summand in the sum C\ 2 +C 22 in the low momentum regime, no(k) ^> 1. This is based on the 
rules of power counting for the degrees of homogeneity h of the linearized collision integral 



which have been given in Ref. 12J. These authors have shown that when the dispersion law 
has the form oj{k) oc k s , and the scattering amplitudes of C\ 2 and C22 scale according to 

M{ek,£k u ek 2 ) = e v M(k, fe 1? k 2 ), (108) 

M(ek,eki,ek 2 ,ek 3 ) = e K M(k,k 1 ,k 2 ,k 3 ), (109) 

the index h is given by 

h[C 12 ] = 5-2r]-d + u, (110) 

h[C 22 ] = 5 - 2k - 2d + 2v, (111) 

where d is the spatial dimension, and v is the exponent in a solution no(k) oc k~ v of C[uq\ = 0. 

Near the critical temperature, the scattering amplitude involving four excitations comes 
from the term g(ip^if)) 2 in the effective Hamiltonian, where ijj vP) creates (destroys) ex- 
citations, while terms as gy/n^ip^ipip give rise to the scattering amplitude involving three 
excitations. In both cases k = rj = since Ai^s oc g,gy/n^. The insertion into Eq. fillip of 
8 — 2, k = 0, d = 3,u = 2 produces h[C 22 ] = 0, while /i[Ci2] = 1. Thus, we can neglect the 
contribution of C22 in the linearized description at low momentum. 

At very low temperature, when the excitations are phonons, the effective description 
is made in terms of a Hamiltonian H[n(x),6(x)] depending on the particle density n and 
the Goldstone mode or phonon field 9 conjugate to n. Now, since each 9- field carries the 
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FIG. 3: U (k,k')k' 3 in the quantum regime, fick 3> 1. The continuous line corresponds to the 
exact kernel for fick = 10. The dashed line is the kernel obtained with Uo(k,k') approximated by 
Eq. (fTT2D for (3ck = 10. 

oscillator factor 1/ ^oo(k), the conjugate field n goes with \Ju){k). Thus, a term in the 
Hamiltonian of order n N produces a scattering amplitude M.^ proportional to Ylf Vk~i, 
where hi denotes the momentum of the quasiparticle i in the process. This means that 
i] = 3/2, according to Eq. ([7]), and k = 4/2 = 2. Therefore, with S = l,u = 1, we find 
MC12] = —4 and /i[C 22 ] = —7, which shows that C\2 prevails over C 22 . 



V. EVOLUTION OF PERTURBATIONS IN THE QUANTUM REGIME f3ck > 1 

In this section we consider the kinetic equation and its solutions when the term Uo(k, k') 
of the collision integral (j22p is approximated in the quantum regime k-^T C cl; <C gn. This 
proceeds by replacing n (l + n ) — > e _/3cfe on the right side of Eq. (12T|) . while Uo(k,k') is 
symmetrically approximated by its leading behavior as fick — > 00, (see Fig. [3]) 

+ e-P ck ' {k - k') 2 6{k' -fc)) , /3cfc^oo. (112) 

The relative size of the second term is 0(e~ l3c( - k '~ k ^). The corresponding subleading contri- 
bution to the collision term may be evaluated by using the Watson's lemma 14 

9 f°° 

' e-P ck {k-k') 2 k' 3 {A(k',t)-A{k,t))dk'~ 
07p~l 3ck k 3 

rdkAfct), /3^0C. (113) 



lQnmn 



8nmn(5 4 c 4 

Although the approximation that arises by ignoring this contribution is necessarily non- 
energy conserving, it may be interesting to explore its consequences, since we would expect 
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this drawback to have a negligible effect as T — > 0. Therefore we consider the evolution 
equation 

d t A(k,t) = — \ f (k-k') 2 k' 3 (A(k',t)-A(k,t))dk' 
lbirmn k J 

9k 5 f k fk'\ 6 (k X <M_ 
~ 16^ J \k) \fr J ^ ' k' 

— A(k,t), (114) 

where, for convenience, we have written the gain term as a convolution. Unsurprisingly we 
see that the loss term reduces to — TB(k)A(k, t) where T^k) = 3/c 5 /(3207rmn) is precisely 
the decay width of Beliaev damping at T = [lOj . 

Although in this regime, n ~ e~^ ck <C 1, we are not in the vicinity of a power law 
spectrum as those appearing in the wave turbulence framework, the above linear opera- 
tor is homogeneous of degree h = —5. Therefore, we can analyze the evolution equation 
with methods similar to the ones we have described before. We define the Mellin function 
corresponding to Eq. fll 14[) 



W(s) = -1 + 60 / v ^ ' x^dx 



x 



s(s 2 - 15s + 74) , . 

Res<4. (115) 



(s-4)(s-5)(s-6) 

The interval of zero winding within the strip of analyticity is (— oo, 0). Rewritten in terms of 
the Mellin-Laplace image J-{s, t) of A(k, t) and the Mellin image *&(s) of the initial condition 
A(k, 0), the evolution equation reads 

AJ^s - 5, A) = v B W{s)F(s, A) + *(s - 5), Re s < 0, (116) 

where v-q = 3/(3207rmn). Now we seek an appropriate solution of the difference equation 
B(s — 5) = —W(s)B(s) in the strip of zero winding, Res G (— oo,0), and we write the 
solution J-(s, A) in the factorized form ( 150]) . To establish this particular solution we require 
that B(s) has neither zeros nor poles for Res G (— oo,0). This requirement is consistent 
with the location of the zero of the Mellin function at s = 0. It must produce only a simple 
pole of B(s) at s = contained in the strip Re G (— oo,0]. The base function turns out to 
be 



r(-f)r(i-z# 



£ \ r / i§ -I- 6 xlk — £ 
s / 1 I io "T L io 5 



m r(^)r(^)r(V) ' (U7) 

and the inverse images of B and 1/B are expressed as Mellin-Barnes integrals along a vertical 
path with Res < 0. The inverse of B(s) may be computed by making the splitting (B(s) — 
1) + 1 which produces a Meijer G function and a delta function: 



Z(*) = 5G$S 



k 



' 10 ' 10 

- — - 

U ' 5' 5 



S{k-l) = Z 1 {k) + 5{k-l). (118) 
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It remains to find the second factor Q(k) of the convolution which gives the solution of the 
initial value problem 

dq 



A(k,t) 



' q 



exp (-v B q 5 t) Q(q)- 



q 



(119) 



In particular, if the initial perturbation is A(k, 0) = C 5(k — ko), the inversion of ^/(s)/B(s) 
yields Q(k) in the form 



Q(k) 



0,3 



7,5 

A 



I _I 

U ' 5' 5 

1 -5+iVn -5-iVn 
' 10 ' 10 



+ C5{k-k ) =Q 1 (k) + C5(k-k ). (120) 



Since all singularities of B(s) and 1/B(s) are located in the half-plane Res > 0, the 
Meijer function G%\(z) vanishes if < z < 1. Some explicit limiting forms are 



^3,3 



1 -1 



G 



0,3 
3,3 



5' 5 
1 -5+iV71 -5-iV71 
' 10 ' 10 

' 10 ' 10 
I -I 

u ' 5' 5 



-M(z-i), z->r 



5cosh(V71vr/10) 



Z —7" OO, 



121} 



122) 



where M ~ 0.48. Consequently, the solution may be rewritten as 



A{k,t) = J q Z i (^j ex P (-VBq 5 t) Q(g)y + exp (-v B kH) Q{k) 

= —Z x ( ^ J exp (-^b&o*) ^(fc - feo) + exp (-^ 5 t) 
&o V fc o / 

+ ^ Z x exp (-^ B g 5 t) Qi(g)y. 



(123) 



It is interesting to find the leading behavior of the last integral as t — > 00. To apply the 
Watson's lemma, we use Eq. (11211) and thus we obtain 



5CM (k . , , , 
zA-)9{k- k ) 



ko 



'ko 

CM 
k 



k 



exp (-v B qH) 



dq 

q 



k \ exp (— Vskyt) 



kr 



vlklH 2 



x9(k-k ), t^oo. (124) 

We see that this integral corresponds to a subleading term of O((T B (k )t)~ 2 ) in the expression 
for A(k,t). Noting that when k > k the second term in Eq. (11231) proportional to Q\{k) 
is subleading with respect to the first one, we obtain the final result for the long time 
asymptotics of the solution 

C f k\ 

A(k,t) ~ —Zx I —J exp (-v B k 5 t) 6(k - k ) 

+C exp (-v B kot)S(k-k ), t-too, (125) 
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or 



A(k,t) 



6C\/lO -2^5 



exp (— fB^o^) ) ^ — > oo, k ^> k . 



(126) 



k cosh(V7l7r/10) 



VI. 



CONCLUDING REMARKS 



In this paper we have studied the temporal evolution of a perturbation of the equilibrium 
distribution of noncondensate atoms in Bose g described by a kinetic equation in- 

cluding only collisions between condensate and noncondensate atoms. Due to the difficulty 
to address this problem without any simplification, we have considered some approxima- 
tions to the kinetic equation in different temperature regimes. The evolution equations that 
arise by approximating the full kinetic equation when the energy of an excitation is small 
compared to the thermal energy turn out to have definite homogeneity properties, and may 
be regarded as kinetic equations for waves in the limit of large occupation numbers. This 
occurs near the critical temperature when the dispersion law is quadratic, and near zero 
temperature for Bogoliubov excitations of very low energy. Notably, by using the explicit 
expression that determines the degree of homogeneity, one may easily show the dominance 
of C\2 over the binary collision term C22 at low momentum. This provides a justification to 
our approach (i. e. dealing only with C12). In the opposite quantum regime where the oc- 
cupation number of noncondensate atoms decreases exponentially with the momentum, we 
recover unexpectedly a kinetic equation of the same type of the wave-turbulence equations 
in the previous situations. 

To obtain the solution of the initial val ue p roblem in all cases, we have applied the 
method developed by Balk and Zakharov 0, Il2| to analyze the behavior of weak-turbulent 
media near Kolmogorov spectra. This is a scarcely used technique due to the difficulties in 
the analytical calculation of the Mellin function and the base function, mainly for binary 
processes 2 -h- 2 in turbulent media. Here, however, this task is made easier as we deal 
with splitting processes 1 2 and the background distribution corresponds to thermal 
equilibrium. 

Just below the critical temperature when the condensate is small, the isotropic pertur- 
bations and the time derivative of the condensate density vanish for long time according to 
power laws, while anisotropic perturbations are unstable. In the thermal regime at low tem- 
perature, the requirement of energy conservation leads to a modification of the naive form 
of the evolution equation, and indeed an additional source term is needed. All solutions 
of the correct equation decay exponentially for large time and the time derivative of the 
condensate density tends to zero according to a power law. In the quantum regime at very 
low temperature the large time behavior of perturbations for the non-condensed dynamics 
turns out to be exponentially decaying due to Beliaev damping processes. 

Finally, one may contemplate the possibility of using a similar analysis to that we have 
presented in this paper to evaluate estimations of the transport coefficients. 
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Appendix: Low momentum analysis of the linearized Boltzmann equation 

In this appendix we give some details about the derivation of the limiting form of the col- 
lision term for small momentum. From Eq. ( fl5|) in the regime near the critical temperature, 
and setting u = cos the following expression for U(k,k') follows easily 

U(k,k') = [n (u(k))[l+n (u(k'))}[l + n (u(k)-u(k'))] 



lQn r a 2 m 2 



m6(k-k')J k'\ n 

-n {u{k) + u(k'))[l + n (u(k))][l + n {u{k'))) 
x^6(u), (A.l) 



where no is the Bose distribution function with zero chemical potential 



"oO") = ^TT- (A.2) 
The coefficients Ui(k, k') are easily evaluated using the standard formula 

Ui(k,k') = — ^— / U(k,k')Pi(u)du. (A.3) 
2 J-i 

To extract the low momentum behavior of Jo{k, q) and Jo{k, q) we first perform exactly 
the integrations of Eqs. ( f28|) and ( f29i) . These may be expressed in closed form in terms 
of polylogarithmic functions Li„(x), where x is an exponential of the some combination of 
reduced energies. As an example, the expressions for J^^i^k, q) near the critical temperature 
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read 



J<(k,q) 



47rn c a 2 



(3m 2 (e^W - l) 2 k 



x 



(3 2 q A -A/3mq 2 In 



1 _ e -j9(w(fc)-w(«)) 



1 _ e -j8(w(fc)+w(g)) 

-8m 2 Li 2 ( e -««(*)-«fo))) - 8m 2 Li 2 ( e -^( fe )+"(*))) 
+ 16m 2 Li 2 (e-^ (fc) )] 6{k-q), 



(A.4) 



J >(M) 



57rn r a 



/3m 2 (e/M fc )-l) 2 £; 



x 



fl _ e -^(«(fl)-w(*))'\ h _ e -/9(u(fc)4«(9)) , \ 

/3 V - 2/3mg 2 In - 



(1 _ e-M*))" 
+4m 2 Li 2 ( e -^(-W-W)) +4 m 2 Li 2 ( e -/%«+"(s))) 
-8m 2 Li 2 (e-^ (9) )]^(g-A;). 



(A.5) 



The form of other kernels J^(k,q) with I ^ are similar. If we substitute the asymptotic 
expansion of Li n (x) near x — 1, 



2 00 /i 

Li 2 (x)~^ + ln(l-x)^^ 



n=i 



n=l 



(1 -a;) n 



x ->■ 1, 



rr 



(A.6) 



we obtain 



647rn r ma 2 1 



/3 2 



A- 3 



In 1 



4 / ' 



64nn c ma 2 1 g 2 + A; 2 
Jo{k,q) ~ 6(q-k) — In 



fe, g -> 0. 



(A.7) 
(A.8) 



/3 2 fc 3 g 2 -fc 2 ' 

These results yield the expression of Hq{x) given in Eq. ( 138]) . once the factor 647m c ma 2 / /3 2 x 
y/ftjm x /3/2m is extracted to define the time scale r in Eq. ( 1361) . 
Near zero temperature the expression for k') reads 



U(k,k') 



9(fc - fc / ) 2 6'(fe - A/; 
327r 2 mn 
+ (A;^A; , )]5(m-1 + 0) 
9(/c + A; n2 



n (o;)[l + n (u')][l + n (w - a;')] 



327r 2 mn 



n (^ W)[l + n (w)][4 + n (u/)]<Kw - 1 + 0). (A.9) 



We note that the 8(u — 1 + 0) arise because the dispersion law u = ck + ak 3 is not strictly 
linear and, since a > 0, the angle between the wave vectors is close to zero. The above 
expression leads directly to the limiting forms given in Eqs. ( 1721) and (1112p . 
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